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ABSTRACT 

Broad-band spectra of the FSRQ (flat-spectrum-radio quasars) detected in 
the high energy gamma-ray band imply that there may be two types of such 
objects: those with steep gamma-ray spectra, hereafter called MeV-blazars, and 
those with flat gamma-ray spectra, GeV-blazars. We demonstrate that this dif- 
ference can be explained in the context of the ERC (external-radiation-Compton) 
model using the same electron injection function. A satisfactory unification is 
reachable, provided that: (a) spectra of GeV-blazars are produced by internal 
shocks formed at the distances where cooling of relativistic electrons in a jet is 
dominated by Comptonization of broad emission lines, whereas spectra of MeV- 
blazars are produced at the distances where cooling of relativistic electrons is 
dominated by Comptonization of near-IR radiation from hot dust; (b) electrons 
are accelerated via a two step process and their injection function takes the form 
of a double power-law, with the break corresponding to the threshold energy for 
the diffusive shock acceleration. Direct predictions of our model are that, on av- 
erage, variability time scales of the MeV-blazars should be longer than variability 
time scales of the GeV-blazars, and that both types of the blazar phenomenon 
can appear in the same object. 

Subject headings: galaxies: quasars: general — galaxies: jets — radiation 
mechanisms: nonthermal — gamma rays: theory — X-rays: general 

1. INTRODUCTION 

The data obtained from the Compton Gamma-ray Observatory (CGRO) mission 
suggest that blazars - the subclass of AGNs which includes FSRQ and BL Lac objects - 
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are strong 7-ray emitters (von Montigny et al. 1995; Mukherjee et al. 1997). In these 
objects 7-ray radiation forms a distinctive spectral component, clearly separated from 
another, lower energy component, presumably produced by the synchrotron process. As it 
was demonstrated by Fossati et al. (1998), blazar spectra form a sequence which can be 
parametrized by their total luminosity. In this sequence, FSRQs are the most luminous 
objects. Both their low and high energy spectral components appear to be the least 
extended to their respective high energies, and their 7-ray luminosity during flares strongly 
dominate over the luminosity of the synchrotron component. The least luminous blazars 
are represented by the X-ray selected BL Lac objects, often called HBLs (High-energy 
peaked BL Lacs). Their low energy spectral component extends up to hard X-rays, and 
their 7-ray spectrum reaches TeV energies. 7-ray flux detected in the TeV-emitting BL 
Lac objects usually does not dominate over synchrotron flux. 

As predicted already in 1978 by Blandford and Rees and later supported by many 
independent observations and theoretical analyses, blazar radiation is most likely produced 
by nonthermal plasma in relativistic jets and Doppler boosted in our direction. Due to 
this, the Doppler enhancement of the jet renders the thermal components of the nucleus 
such as the UV radiation from an accretion disc, the X-ray radiation from the disc corona, 
and the infrared radiation from the hot dust - all presumably emitted isotropically - 
significantly diluted. But, at least in FSRQs, the optical- UV broad emission lines (BEL) 
are clearly detectable. They provide information about the redshift as well as the energetics 
of the central engine and about radiative environment of the sub-par sec/parsec scale 
jets. It turns out that the energy density of the BEL region on sub-par sec/parsec scales 
as measured in the comoving frame of the jet usually dominates the energy density of 
the synchrotron radiation. Within the framework of the internal shock model for the 
generation of radiative outbursts, that implies that the 1-10 day time scale 7-ray flares 
are produced by Comptonization of emission lines (external-radiation Compton [ERC] 
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model) (Sikora, Begelman, & Rees 1994), rather than by Comptonization of synchrotron 
radiation (synchrotron-self-Compton [SSC] model), the latter presumably dominating the 
7-ray production in the low luminosity BL Lac objects. 

The ERC model successfully predicts the observed location of the spectral break in 
the high energy spectra of FSRQs to be in the 1-30 MeV range without the necessity of 
postulating a break in the power-law electron injection function. According to the model, 
this spectral break simply reflects the break in the electron energy distribution caused by 
the effect of particle cooling (Blazejowski et al. 2000; Sikora et al. 2001; see also §2). If this 
interpretation is correct, the change of the spectral slope around the break resulting from 
cooling should be unique and equal to Aa XJ = 0.5, where a : F u oc v~ a . In most FSRQs, 
Aa xl < 0.5, and the values of Aa xl lower than 0.5 can be explained by invoking the 
possible dominant contribution of the SSC process in the soft /mid X-ray band. However, 
the so-called MeV-blazars present a challenge for the model. Their X-ray spectra are very 
hard, a x = 0.3 — 0.5, and the high-energy 7-ray spectra are very soft, 07 > 1.4, and, 
therefore, Aa xl > 1 (Blom et al. 1998; Bloemen et al. 1995; Tavecchio et al. 2000). There 
were several attempts to explain such spectra in terms of the ERC model. Sikora et al. 
(1997) proposed an inhomogeneous model, according to which the X-ray spectra are due 
to superposition of multiple ERC components produced at different distances in a jet and 
having low energy spectral cutoffs at frequencies determined by the plasma temperature, 
assumed to increase with a distance. This model, however, cannot reproduce the correlation 
of the X-ray and 7-ray light curves during the 1996 February flare observed in 3C 279 
(Wehrle et al. 1998; Lawson et al. 1999). Georganopoulos, Kirk, & Mastichiadis (2001) 
suggested that the steep 7-ray spectra can result from the fact that they are softened by 
the Klein-Nishina effect. But this effect is important in the EGRET band (30 MeV - 3 
GeV) only if energy of the seed photons involved in the ERC process is more than 10 times 
greater than the typical energy of the BEL photons. Finally, Blazejowski et al. (2000) 
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demonstrated that soft 7-ray spectra can result from superposition of two components, 
one produced by Comptonization of BEL, and another due to Comptonization of near-IR 
radiation produced by hot dust. This model requires the high energy cutoff of the ERC(IR) 
component to be in the range 3-30 MeV and such a condition can be satisfied only via 
fine-tuning of the maximum electron energy to be within the range 3 x 10 3 — 10 4 m e c 2 . 
Another weakness of this model is that the spectra superposed from the two components 
do not provide a good fit of the approximately power-law spectra observed in the EGRET 
band. 

In the present work we assume that electrons are accelerated via a two-step process 
and that their injection function takes the form of a double power-law with a break at the 
energy which divides the regimes of dominance of two different acceleration mechanisms. 
Here the division energy is the threshold energy for the resonant scattering of electrons by 
Alfven waves (required to have efficient diffusive shock acceleration of electrons). With this, 
we demonstrate that the MeV-blazar type spectra are produced when 7-rays result from 
Comptonization of the infrared dust radiation, whereas the GeV-blazar type spectra are 
produced if the 7-ray flux is dominated by Comptonization of BEL. Our work is organized 
as follows: in §2, we discuss the cooling effect and demonstrate that for a distance range 
corresponding to the observed variability time scales, the cooling break is located within 
the 1-30 MeV range, as implied by observations. In §3, we discuss the issue of the electron 
acceleration efficiency and present the motivation for introducing the double power-law 
approximation for the electron injection function. Such an approximation is used in §4 
to calculate the time-averaged spectra produced at different distances in a jet by shocks. 
There, we also demonstrate that adopting the same electron injection function, one can 
reproduce the typical spectra of both MeV- and GeV-blazars. Our results are summarized 
in §5. 
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2. ELECTROMAGNETIC SPECTRA 

2.1. Cooling effect 

The basic feature of the high energy spectra in FSRQs — a spectral break between the 
X-ray and the 7-ray bands — has a natural explanation in terms of the ERC model. In 
this model, X-ray spectra are produced by electrons with radiative cooling time scale t' cool , 
longer than the source life-time t'^ (slow cooling regime), whereas 7-rays are produced by 
electrons with t' cool < t'^ (fast cooling regime). Noting that the angle- averaged cooling rate 
of electrons, when dominated by Comptonization of external radiation, is 

\l\ - — »W7 , (1) 

where 7 = E' el /m e c 2 is the random Lorentz factor of the electron and u' ext is the energy 
density of the diffuse external radiation field, both as measured in the source comoving 
frame, we obtain the angle-averaged cooling time scale of the electron to be 



, 7 m e c 1 , . 

^cooi — rn" — / 5 (.2 j 



where u' ext is the energy density of an external radiation field. Then, from t' cool = t'^ = Vtfi, 
where tfi is the observed time scale of the flare, and 

T> = \ (r>) 

T(l-f3cos6 obs ) 1 ; 

is the bulk Lorentz factor, the break in the electron distribution is located at the energy 

m e c 1 

T U ext tfiV 

For 7 < 7c, the slope of the electron distribution is the same as the slope of the injection 
function; for 7 > j c , the slope of the electron energy distribution is steeper by As = 1 
(s : iV 7 oc 7" s ). 



-7- 



In the internal shock model, the lifetime of the flare is equal to the lifetime of the 
shock(s) and that is equal to the collision time scale of the two inhomogeneities, which prior 
to the collision are assumed to propagate down the jet with different velocities. Thus, we 
have the distance range where the shock is active to be 

A rcoU = ctflVF , (5) 

and because time scales of flares in FSRQs are rarely shorter than 1 day, those flares are 
most likely produced at distances larger than 0.1 parsec. At such distances, the contribution 
to u' ext is dominated by BEL and infrared radiation from hot dust. Noting that 



U ext 



- c J I' ext dQ' = -J I ext V- n 2 dQ ~ u dtff T 2 (6) 



where 



Vm r(l-/3cosU (7) 
and 9 in is the angle between the photon direction and the jet axis, we predict that the break 

in an electron energy distribution at 7 C should correspond to a break in the electromagnetic 

spectrum at frequency 

, 2 , _ (m e c\ 2 is diff ( m e c 2 \ 2 is diff (V x z 



\°t J u ext t fl \ a T J u diff /\r coU \Y J 
and that the spectrum should change the slope around u c by Aa xl ~ 0.5 (<— As = 1). 



2.2. External radiation fields 

Broad emission line region 

According to the interpretation of the data obtained in many reverberation campaigns, 
production of broad emission lines in quasars is stratified and peaked around a distance 
(see, e.g., Peterson 1993; Kaspi 2000; Sulentic, Marziani & Dultzin-Hacyan 2000) 

tbel ~ 3.0 x 10 17 v/Z^cm (9) 
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Unfortunately, the detailed dependence of line luminosities on distance and geometry of the 
BEL region is poorly known. For our illustrative purposes we assume that the BEL region 
is spherical and that the fraction of the central UV radiation reprocessed into lines at a 
distance r is £,bel{t > r BEL ) = d£ BEL /d\nr oc (r BEL /r) q and £bel{t < r BEL ) = 0. With 
that, the energy density of BEL field, as measured in the co moving frame of the radiating 
plasma, can be approximated by 



where L BEL = £ BEL L UV and £ BEL = j™ EL (d£ BEL /dr)dr. Note that for r > r BEL , the 
value of ubel(j') is dominated by radiation coming from smaller distances, from around 
tbel- However, since this radiation is redshifted in the source frame, the contribution to 
u' BEL (r) is dominated by the fraction of ubel(j') which is determined by broad emission 
lines produced at a distance r. 

Hot dust 

Evidence that hot dust is present in blazars is indirect. Infrared emission of hot dust 
is directly measured in quasars which are observed at larger angles to the jet axis than 
blazars. The ratio of the IR flux to the UV flux in these objects shows that the fraction of 
the UV radiation reprocessed by the dust into the infrared, i.e. the dust covering factor 
£m, is of the same order as £ B el (Sanders et al. 1989). Dust is probably concentrated in 
molecular tori, but optical extinction suggests that its distribution in the normal direction 
to the equatorial plane doesn't have any sharp boundary (Baker 1997). One of the unknown 
aspects of dust in quasars is the minimal distance from the nucleus where it can exist. This 
distance can be limited by the maximum temperature that dust can survive, 



where T d>max ~ 1500K, or by the inner edge of the torus if it is larger than r d (T d max ) (see, 
e.g., Yi, Field, & Blackman 1994). As observations of individual objects show, there is a 




(10) 




(11) 
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very large scatter in the dust amount and its distance distribution amongst various quasars 
(Polletta et al. 2000; Andreani et al. 2002). 

In the present calculations we approximate the dust distribution by assuming that it is 
spherical and enclosed within a given distance range with a constant covering factor. With 
these assumptions the energy density due to dust, as measured in the source comoving 
frame, can be estimated using the formula 



where A = ln(r dtTnax /r djmin ) ~ 5. The dependence of u' uv and u' IR on r is illustrated in Fig. 
1. For such radiation fields we plot the dependence of 7 C on r on Fig. 2 and of v c on r Fig. 
3. It is apparent from the Fig. 3 that for a very wide range of r, the break hu c is located 
within the photon energy range 1 MeV - 30 MeV, in agreement with observations. 



As it was demonstrated by Sikora & Madejski (2000), the energy flux in powerful jets 
in quasars cannot be dominated by pair plasma. This is because such a jet would produce 
much larger flux of soft X-ray radiation than is observed. Hence, we assume that the inertia 
of the jet is dominated by protons. In this case, the structure of shocks and the structure 
of Alfven waves generated around the shocks are both determined by protons. Being 
resonantly scattered by these waves, the protons jump back and forth across the shock front 
and participate in the 1st order Fermi acceleration process (Bell 1978; Blandford & Ostriker 
1978). Those protons which reach energies > 10 9 GeV interact efficiently with ambient 
photons and trigger (mainly via photo-meson process) synchrotron-supported pair cascades 
(Mannheim & Biermann 1992). However, such a model fails to reproduce the very hard 
X-ray spectra of FSRQs (Sikora & Madejski 2000). This may indicate that in the context 




1 + (r/r d:Tnin ) 2 A ' 



1 1 



(12) 



2.3. Electron injection function 



-10- 



of those models, too few protons reach sufficiently high energies to power pair cascades. 

An alternative possibility is that the high energy radiation in FSRQs is produced 
by directly accelerated electrons. However, for efficient acceleration of electrons by 
the Fermi process via resonant scattering off Alfven waves, the electrons must be first 
preheated/preaccelerated up to energies 7^, at which point the magnetic rigidity of electrons 
becomes comparable with magnetic rigidity of thermal protons, i.e. when their momenta 
are equal: 



is the average thermal proton energy in the shocked plasma, k is the amount of energy 
dissipated per proton in units of m p c 2 , and rj Ptt h is the fraction of the dissipated energy 
tapped to heat the protons. In the case of the two intrinsically identical inhomogeneities 
(Sikora & Madejski 2001) 

2(ry ri )V2 • (15 

where > Ti ^> 1 are bulk Lorentz factors of inhomogeneities prior to the collision. For 
the reasonable assumption that r 2 /Ti < 10, the thermal proton plasma is at most mildly 
relativistic, i.e. 7 Pjt / l — 1 < 1. 

Noting that 

n e m e {^ inj - 1) ~ n p m p r] e K = n v m v -^—{^ t h - 1) , (16 

Vp,th 

where 7j n j is the average Lorentz factor of injected electrons and i] e is the fraction of the 
dissipated energy used to accelerate electrons, we find that for 7^ and j in j 1 




(13) 



where 



l P ,th - 1 = rj P ,thK 



(14) 




(17) 
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Hence, for non- or mildly relativistic shocks and for r} Pjt h ~ Ve, the threshold energy for the 
diffusive shock acceleration of electrons significantly exceeds the average energy of electrons, 
even if n e = n p . This implies that the often-considered bulk preheating process is not able 
to provide an adequate number of electrons with 7 > jp. However, as it was recently 
demonstrated in numerical PIC (particle-in-cell) simulations, the preheating/ preacceleration 
mechanism can have a stochastic character and a non- negligible fraction of electrons can 
reach such energies (Dieckmann et al. 2000; Shimada & Hoshino 2000). Furthermore, the 
fact that even in blazars with the hardest X-ray spectra, the X-ray spectral indices ax 
are always greater than indicates that the largest number of electrons is injected at low 
energies, and, therefore, that there is no evidence for a bulk preheating process forcing most 
of electrons into equipartition with protons. This strongly suggests that in a similarity to 
the diffusive shock acceleration operating at 7 > 7^, the preacceleration process also has 
stochastic character and injects electrons with a power-law energy distribution. Motivated 
by this, we assume that electrons, being accelerated by the two-step process, are injected 
with the double power-law energy distribution 



where Ch = Cffp Pl . The injection break at 7^ given by Eqs. (13) - (15) results in a break 
in the electromagnetic spectrum at the frequency 




(18) 



(19) 



3. 



MeV-BLAZARS vs. GeV-BLAZARS 



For 2 < r 2 /r! < 10 and 1/3 < r] Pjt h < 1/2, the break in the electron injection function 
is 7i? ~ 10 3 (see Eqs. 13 - 15 and Fig. 2). Another break appears in the electron energy 
distribution at 7 C (Eq. 4). If the cooling of electrons is dominated by Comptonization 
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of BEL, then 7 C <C 7^, and the frequency range [v c , v B \ overlaps significantly with the 
EGRET band. Hence, the spectrum in the EGRET band should show a slope a ~ pi/2 and 
should be steeper there than in the X-ray band by ~ 0.5 (or less, if the contribution of the 
SSC radiation into X-ray band is taken into account). For pi < 2.0 such spectra become 
representative for short term flares in GeV-blazars (Pohl et al. 1997). 

At larger distances, the production of 7-rays is dominated by Comptonization of 
infrared radiation of hot dust. Luminosity of that radiation is comparable or even larger 
than the luminosity of the broad emission line light, but because of greater distance, its 
energy density is much smaller than energy density of the diffuse radiation in the BEL 
region. As a consequence, the cooling break energy, 7 C , is now much larger, approaching 
the value comparable to the value of 7^ (see Fig. 2). Hence, there is no longer an extended 
spectral plateau with a — pi/2, and both u c and v F conspire to produce the spectral break 
Aa X7 = 0.5 + {ph — Pi)/2. Since now the energy of the Comptonized diffuse photons is 
virIvbel times lower than in the BEL region, an approximate location of this break is at 
a frequency which is by that factor lower than vp in the BEL region. Hence, the spectra 
produced in the EGRET band at larger distances are predicted to have a spectral slope 
~ a 3 = Ph/2. For pi < 2.0 and p h > 2.4 such spectra become representative for MeV-blazars 
(Tavecchio et al. 2000). 

The above scheme can also explain two other observed differences between MeV-blazars 
and GeV-blazars. One is that the X-ray spectra are generally much harder in MeV-blazars. 
In the context of the scenario above, this presumably results from the lower contribution of 
the SSC component to the X-ray band in the MeV-blazars than in the GeV-blazars, which 
is the case, e.g., if (u' IR /u' B ) r2 > (u' BEL /u' B ) r2 , where u' B is the magnetic energy density and 
r 2 3> 77 ~ r BB L- The other difference is that in the spectra of MeV-blazars, in contrast 
to GeV-blazars, the thermal UV bumps are quite prominent (Tavecchio et al. 2000). This 
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difference can be explained by noting that magnetic fields are weaker at larger distances, 
and therefore in MeV-blazars, the synchrotron spectra are shifted to lower frequencies, 
revealing the UV bump. All of the above effects are illustrated in Figs. 4 & 5, where 
we present the time averaged spectra of radiation produced at two different distances by 
electrons injected with the same energy distribution, Q, given by Eq. (18). The spectra are 
computed using the internal shock model where the shock propagates down the conical jet, 
and including the following dominant radiation processes: synchrotron radiation - SYN; 
Comptonization of synchrotron radiation - SSC; Comptonization of BELs - ERC(BEL); 
and Comptonization of infrared radiation of hot dust - ERC(IR). The model is presented 
in Blazejowski et al. (2000) and Sikora et al. (2001); for a more comprehensive description 
of the model see Moderski et al., in preparation. The model input parameters used to 
calculate the presented spectra are specified in the figure captions. 

In order to illustrate the relative location of the breaks at 7 C and 7^ in these two 
specific models, we present in Figs. 6 and 7 the energy distributions of electrons. They are 
computed at a distance of the shock termination (which in our models is at r = 2ro) and 
presented together with the electron injection function to demonstrate the cooling effect. 

It should be emphasized here that the value of 7^ does not depend on how the fraction 
(1 — rj p ,th) of the energy dissipated in the shock is shared by other consistuents than thermal 
protons (see Eqs. 13 - 15). We also note that the comparison of the FSRQ spectra during 
outbursts against the spectra calculated from our model indicates that the best matching 
is provided by the models where the electron energy density exceeds the magnetic energy 
density by a factor ranging from a few up to more than 10. 
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4. CONCLUSIONS 

We demonstrated that the appearance of a FSRQ as a GeV-blazar or a MeV-blazar 
can be explained as depending on the distance of the maximum rate of energy dissipation 
in a jet. The GeV-flat spectra are predicted by our scenario to originate at distances 
where the production of 7-rays is dominated by Comptonization of BEL, while the spectra 
characteristic of MeV-blazars are predicted to originate at distances where the production 
of 7-rays is dominated by Comptonization of hot dust radiation. An additional ingredient, 
required to reach satisfactory unification between these objects, is the assumption that 
electrons are injected with the double power-law energy distribution, with the break at 
7 ~ 10 3 . Such location of the break coincides numerically with the threshold energy for 
the diffusive shock acceleration of electrons. Below that energy, the electrons must be 
accelerated by a different mechanism, e.g. following instabilities driven by shock-reflected 
ions (Hoshino et al. 1992; McClements et al. 1997; Shimada & Hoshino 2000; Dieckmann 
et al. 2000); or following reconnection of magnetic fields (Romanova & Lovelace 1992; 
Blackman 1996). 

As observations of X-ray spectra in MeV-blazars show, less energetic relativistic 
electrons are accelerated with the power-law energy distribution, with index pi < 2. Future 
studies of dependence of the spectral shape on time scales of outbursts can be used to 
verify our idea that spectral breaks Aa xl > 1 result from superposition of two breaks, one 
produced by the cooling effect and another one reflecting the break in the electron injection 
function at the threshold energy for the Fermi acceleration process. 

It is very likely that some of the MeV-blazar phenomena can be produced at distances 
where the source is already optically thin at mm-wavelengths. If this is the case, the model 
predicts a correlation between variability of the flux around the synchrotron-self-absorption 
break at mm-wavelengths and in the 7-ray band, on typical time scales of a month. This 
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possibility is supported by recent polarization measurements at very high radio frequencies. 
Those observations reveal that magnetic field in the very compact radio cores is dominated 
by its prependicular component and this is consistent with the model of the radiation 
produced in situ by transverse shocks (Lister 2001). Depending on the amount and 
covering factor of the hot dust in the surrounding medium, such shocks can be sites of the 
MeV-peaked 7-ray blazars. 

Finally, we would like to comment that the MeV-blazar and GeV-blazar phenomena 
can interchangeably appear within the same object, as likely is the case in PKS 0208-512 
(Stacy et al. 1996; Blom et al. 1996). With such objects, correlation between the spectral 
type and the variability time scales can be studied directly, without the bias introduced, for 
instance, by the dependence of the variability time scales on the black hole mass which can 
be very different from object to object. We note that the future sensitive 7-ray telescopes 
such as GLAST will provide us with excellent spectrally resolved light curves and will 
significantly constrain the applicability of radiative and particle acceleration models for 
those enigmatic sources. 

This project was partially supported by Polish KBN grant 5 P03D 002 21 and a 
Chandra grant from NASA to Stanford University via the SAO award no. GOO- 1038 A. 
M.S. thanks SLAC for its hospitality during the completion of part of this work. 
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FIGURE CAPTIONS 

Fig. I. — The dependence of the energy densities of external radiation fields, u' BEL (dashed 
line) and u' IR (solid line), and of energy density of magnetic field, u' B = B' 2 /8it (dotted line) 
on the distance from the central source, r. The curves are calculated using Eqs. (10) and 
(12) for: r BEL = 5.2 x 10 17 cm; r d , min = 2.6 x 10 19 cm; T = 15; L BEL = 3.0 x 10 45 erg s" 1 ; 
L d IR = 9.0 x 10 45 erg s _1 ; q = 1.0; A = 3.2; and B' = (1.2 x 10 18 cm/r) Gauss. 

Fig. 2. — The dependence of the electron energy distribution breaks, 7 C and 7^, on distance. 
The 7 C curves are calculated using Eqs. (4) and (5) for r = Ar co u: ^ C (UV) (dashed line) - 
assuming that external radiation field is totally dominated by u' BEL ; and / ~f c (IR) (solid line) 
- assuming that external radiation field is totally dominated by u' IR . The value of 7^ is 
obtained from Eqs. (13) - (15) for Y 2 /Vi = 3.0 and rf Pjt h = 0.5. All other model parameters 
are the same as in Fig. 1. 

Fig. 3. — The dependence of the "cooling break" , v c on the distance from the central source r. 
The break is calculated for the same model parameters, as in previous figures and assuming 

e obs = i/v (-z> = r). 

Fig. 4. — The model spectrum of the time averaged flare produced by the shock formed 
at a distance r = 6.0 x 10 17 cm and terminated at r = 2r . The shock propagates down 
the conical jet with a half-opening angle 9j = 0.02. The electrons/positrons are injected at 
the rate given by Eq. (18) for Ci = 3 x 10 49 s _1 , pi = 1.8, and ph = 2.8. All other model 
parameters are the same as in previous figures. The two areas confined by vertical dashed 
lines represent the 2-10 keV and 30 MeV - 3 GeV bands, respectively. The sharp feature 
around v ~ T 2 u BEL is the artifact of two approximations of our model: first, that we do not 
take into account upscattering of external photons produced at the distances smaller than 
the actual position of the shock in a jet - which would contribute to spectrum at v < T 2 ubel, 
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and second, that we use the average BEL frequency. 

Fig. 5. — Same as Fig. 4, but for the shock formed at a distance r = 2.0 x 10 19 cm. 

Fig. 6. — Energy distribution of electrons 7 2 A r 7 at r = 2r for r = 6.0 x 10 17 cm, shown as 
a dotted line. For illustration, we also show the electron injection function multiplied by the 
period of the shock operation ^Qt'^ as a solid line. 

Fig. 7. — Same as Fig. 6, but for the shock formation distance r = 2.0 x 10 19 cm. 



-21- 




- 22 - 




-23- 




17 



18 19 

log r [cm] 



20 



-24- 




- 25 - 




-26- 




12 3 4 

log y 



-27- 




12 3 4 

log y 



